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We define a Hermitian phase operator for zero mass spin one particles (photons) by taking account 
polarization. The Hilbert space includes the positive helicity states and negative helicity states with 
opposite circular polarization. We define an operator which corresponds to the physical process of 
reversing the sense of polarization and acts as a bridge between positive helicity states and negative 
helicity states. The exponential phase operator obtained using the entire set is unitary and acts 
as ladder operator over all the states. The phase operator derived from this exponential operator 
satisfies the canonical commutation relations with the number operator. We have calculated the 
density matrix and the phase probability distribution of various states like coherent states, squeezed 
states and thermal states, to illustrate the utility of our operator. 
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I. INTRODUCTION 



As almost all measurements on electromagnetic waves are measurements of intensity or phase; there has been an 
abiding interest in a phase operator for light. Having identified the number operator as the quantity proportional 
to intensity, Dirac [T] decomposed the annihilation operator for the harmonic oscillator as a = E\fti = e l ^\ffi. Here 
4> is the phase operator and h the number operator of the oscillator. The number operator was to be the quantum 



h,<j> 



i. As he pointed 



equivalent of the classical action and satisfy canonical commutation relations with the phase, 

out, this phase operator is not well defined as there is a lower bound on the eigenvalues of the number operator. In 
addition, the phase is periodic and defined modulo 2ir. Subsequent workers [2THT] all had the same problems of either 
lack of unitarity of the exponential operator or failure to obtain the canonical commutation relations for the number 
operator with the phase operator. An excellent review of the above ideas as well as those to be described below was 
given by Lynch |12j . 

Pegg and Barnett [SI U\ suggested the most popular way round these difficulties. The idea is to truncate the 
oscillator states to a finite number say s. The operators in the resulting finite dimensional Hilbert space are well 
behaved, this is extremely useful as a practical tool. However, questions have been raised [5HTU] of whether such 
a weak convergence to the limiting infinite dimensional operator gives the same results as the calculations directly 
using infinite dimensional operators. In view of these difficulties, Noh, Fougeres and Mandel [131 IH] even suggested 
that a unique phase operator should not be sought. They suggest that we should only use an operational definition 
of phase based on experiment performed. 

With the recent development of quantum tomographic methods to study entanglement and deduce the Wigner 
distribution, the phase or more accurately the phase difference has again become very topical [15] . There has been 
many subsequent efforts to define a phase operator both for photons and for the closely related problem of the 
Harmonic oscillator JTBHH] . Many of these, especially the work by Ban, Luis and Hradil, use products of two boson 
operators which are especially well suited to study parametric processes and higher order interference but are also 
used to cover the case of two orthogonal light polarization modes. These two boson representations are best suited to 
treat problems where the total number of photons is constant, as the two boson representations used are isomorphic to 
the angular momentum group. However as pointed out by Bjork et. al., [22] these are not well defined in the vacuum 
state and indeed require a separate definition of the two mode vacuum state, which is annihilated by all elements of 
the Schwinger two boson algebra. Further, the addition of one more photon which can take place spontaneously or 
by stimulated emission, with very different consequences for the phase of the state is not naturally incorporated in 
the theory. 

There is another representation of the electromagnetic field, with the photon helicity as an internal quantum number, 
which was advocated by Bialynicki-Birula |23| and independently by Sipe |24j . this suggestion made in the context 
of the photon wave function was subsequently shown to yield many other properties of the electromagnetic field in a 
fashion such that the passage from the classical theory to the quantum electrodynamics is particularly transperant, 
[2"5H2"5] . Our representation of the field in this paper is also over 2 ® (0,oo), rather than the usual product state 
formalism. 

The Hilbert space we consider is a tensor product of a two internal dimensional space of polarization and the Fock 
space of a single polarization mode. This agrees with the representations found by Wigner, [301 HU f° r massless 
relativistic particles. Bargmann and Wigner showed the appropriate relativistic representation for photons is two 
dimensional, this was elaborated later |32j to show that it is indeed the natural representation to take into account 
all polarization phenomena. As the detailed calculations of Kim et. al. show the operators for the commonly used 
devices like polarizers, phase plates are particularly simple in this representation being constant matrices. The action 
of such operators is similar on states with the same polarization but different number of photons. 

As the Riemann-Silberstein vector E + iB as the photon wavefunction advocated by Bialynicki-Birula and Sipe 
yields a representation for photons in both position and momentum space, this method can be also used to study 
nonlocal correlations of the EPR type, which are not easily defined in the usual product space representations. Also it 
is known that the the Cartesian components of the light polarization need to commute to give a definition consistent 
with Maxwell's equations [33H35] . 

Sperling and Vogel [36] recently pointed out that removal of a right circularly polarized photon is the same as 
the addition of a left circularly polarized photon. They then suggest treating the left circularly polarized photons as 
antiparticles of the right circularly polarized photons and used this to construct a unitary Susskind Glogower operator 
E. The difficulty with this theory is that there is a difference in angular momentum between the vacuum for left 
circularly polarized photons and that for right circularly polarized photons. 

We suggest a modification of the Sperling- Vogel procedure, based on two requirements. 
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First, the quantum phase being an angle, the corresponding momentum must be an action variable. The number 
operator acting on a Fock state then measures the number of action quanta in the state. No experiment directly 
measures the action. However, the number of action quanta is equal to the number of azimuthal angular momenta in 
that state, which is measurable; so we have a suitable proxy. The essence of our suggestion is that action and phase 
form a conjugate pair while energy and time form another conjugate pair and treating them as freely exchangeable 
may be part of the reason for the earlier difficulties. 

Second, by passage through a phase plate, it is possible to transform a state of definite momentum and polarization 
"a mode" , into one of the orthogonal polarization, "an orthogonal mode" , while preserving the phase or changing it 
trivially. So the Hilbert space relevant to phase measurements is that spanned by both the polarization modes taken 
together. 

If the number of right circular polarized quanta present is denoted by n + > then the number of left circular 
quanta present is obviously rt_ < or vice versa. This method of labelling of the Helicity eigenstates is consistent 
with the requirement that if the right circularly polarized photons have positive values of azimuthal polarization and 
positive action, then the left circularly polarized photons have negative values of azimuthal angular momentum and 
consequently negative values of action. Our method of labelling permits a natural extension to cover the negative 
integers required to attain unitarity. Newton |4j, suggested the use of such a two dimensional representation, but in 
terms of energy eigenstates, consequently he had to suppress the negative energy states which were only a mathematical 
intermediate step in his model. The angular momentum of light both orbital and spin has been measured and shown 
to add algebraically [37]. 

In the next two sections, the Hilbert space and the Susskind Glogower operator are defined and the operator is 
shown to be unitary, while the corresponding phase operator is Hermitian. We then find the eigen states of the 
Susskind Glogower operator and show they are orthogonal but non-normalizable. All the previous results cited above 
can be recovered by restricting our space to states of a single polarization. While we outline our results using the two 
opposite circular polarizations, they hold in any polarization basis. Finally, using the projectors on to the different 
polarizations we construct operators which change the number of photons in one polarization state while leaving 
the other unchanged. As this operator for photon removal (addition) commutes with the corresponding annihilation 
(creation) operator acting on the whole space, we can use these to construct any arbitrary photon state. 



II. THE PHASE OPERATOR 



The positive frequency components of a single momentum mode of the electromagnetic field are usually written as 
[38], 

E(r,t) cx a+e 1 ^^ + a^e 1 ^'^. (1) 

Here a + is the annihilation operator for right circularly polarized light, and a_ is the annihilation operator for left 
circularly polarized light. These operators satisfy the commutation relations, 



a + ,a\ 
[a+,a-\ = 



it 



a+, a 



= 1, 
= 0. 



If P is the helicity operator then 



P\n± >= ±|n± > and (P) 2 = 1 



(2) 



(3) 



As there are two states for a given momentum k and these are orthogonal (Malus law) the natural representation 
appears to be in terms of a two-dimensional internal space. 
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We use the operator P to define operators which project into the circular polarization eigenstates. 
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(6) 



To incorporate the polarization directly into our description, we note that the commutation relations for the annihi- 
lation and creation operators imply that if \n > is an eigenstate of the number operator with eigenvalue n then 

a) a(a\n >) = (n — l)(a\n >), 

a)a{a)\n >) = (n + l)(a)\n >). (7) 

Both are eigenstates with eigenvalues n — 1 and n + 1. XTie commutation relations impose no further restrictions on 
n. The series termination requirement yields the fact that n must be an integer. However the set of all non-negative 
integers as well as the set of all non-positive integers both satisfy equation Q. 



a + \n > 
it 



1 >, 



a' + \n > = vn+T |n + 1 >, 



n = 0, +1, +2, +3 . 
//! in - 1 >, 



a+\0 >= 0, 

a_ \n > 

al|n > = y/\n + 1| |n + 1 >, 
o^|-l>=0, =^n = -l,-2,-3.... 
We now apply this result to our formulation; the modified annihilation operator of system is written as 



it 



o 

it 



(8) 



(9) 



(10) 



Where a m + acts on positive states and a m _ acts on negative states and are given by, 



l m-\- 



= P 




= P 



- 1 ><o| 




i + p 



= 0, 



a-,P 



= 0. 



(11) 



For positive helicity states a + is the annihilation operator and for negative helicity states a\_ is the annihilation 
operator and a v is an operator which serves as bridge between two subspaces of positive helicity states and negative 
helicity states. It causes the vacuum of positive helicity states to change into vacuum of negative helicity states and 
vice- versa. Note that quantum theory requires them both to have the same energy at least in the absence of a medium. 
The state space is spanned by the basis, 



\n+ > 




and 



with, 



\n+ >= 



1 > 



/n^\n + - 1 >, a m +\n- >= 0, 
n+ =0,1,2,3,... 



(12) 



(13) 



■,t 



a,; ri _|n_ >= vjnTj|n_ - 1 >, 



it 



% _\n+ >=0, 



1,-2,-3,. 



(14) 



Thus our annihilation and creation operators are 



a+ 
1 >< 01 a_ 



and 



al |0 >< -1| 




al 



(15) 



(16) 



Using these annihilation and creation operators, number operator and an exponential of the phase operator arc 
constructed in a similar way as done by Susskind and Glogowcr. The number operator in this representation is, 



and the exponential of the phase operator 



h+ 
h--l 



(17) 



(18) 



The exponential of phase operator was found to act as ladder operator over all the states. 

E m \n m >= \n m - 1 > , V n m = n + and n_, 
E m \0 >= | - 1 > and ££j - 1 >= |0 > . 

The exponential of phase operator is unitary as, 

E — E =1 



(19) 



(20) 



so we can write E m = e l ^ m and hence the phase operator is Hermitian and satisfies the canonical commutation 
relations with number operator as 



E m so 



'm 5 



(21) 



Therefore, the canonical representation for the number operator in the basis of phase ((f) eigenstates is written as 

. d 



dej)' 



The commutation relation of our modified annihilation and creation operators becomes, 

[ct m , 



P-|0><0| 

o P + l-ix-il 



(22) 



(23) 



So the non-diagonal matrix elements vanish and diagonal matrix elements give +1 in the positive number states and 
-1 in the negative number states as required. Most interestingly they vanish in both vacuum states 



< n + \ |a TO ,aLJ \n+ >= 1, 

< n_| [a m ,«L] \ n - >= - 1 ! 
< -1| [a m , a\ n ]\-l >=< 0| [a m , a) m ] |0 >= 0. 



(24) 



This means all the commutation relations vanish between the vacuum states, showing that there is no energy difference 
between the vacuum states of either polarization. The two vacuum states are degenerate. So the appropriate vacuum 
state for bosons is their symmetric combination as required by quantum theory. The eigenvalue equation for E m gives 
the phase states 
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Thus we conclude that the two vacuum states differ in phase and are orthogonal 
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III. PHASE PROBABILITY DENSITY 



Using this phase state, we next calculated probability of phase distribution in coherent states, squeezed states 
and thermal states. If the density matrix for the pure phase state \ip > is given by p, then the phase probability 
distribution p p (cj)) is given by 

p p (4>) =< ip\p4,\ip > where P<p = \<j> m >< <f> m \ 

-l 
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(27) 



where 4> is the phase and n is any integer between and infinity. The density matrix for phase state is given by 
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The corresponding probability density is 
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(29) 



The first two terms give phase probability distributions of the right circular polarization and left circular polarization 
and the remaining two terms give the polarization interference. 



A. Phase probability distribution for a coherent state 

Taking account the polarization, the coherent state is written as 



a >= e 2 



E 
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a" 





n - 1 > 



(30) 



Positive value of n represents all right circular polarization states and negative value of n represents all left circular 
polarization states. Substituting this value in the expression of phase probability density gives, 
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(31) 
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Solving for all the terms and taking a — 1 we get the phase probability density for right circular and left circular 
polarization as, 



oo oo 



R(M - -Lp-N* >T >T a a r i(n-m) 



(4>) = E E ^e"*—)* (32) 



while the polarization interference term is 

= J_ e -i-i 2 y y 5_^ e »(n+ m )0 + _l e -H 2 y y ELJLe-nr*™)* (33) 

^/Av; 2 Z^ Z^ ^Tf^T 2tt Z^ Z^ Jrf^x v ; 

We plot phase probability density of all the polarization states with the phase in the limit — 7r to +tt taking a = 1, 
We get the plots as, 
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FIG. 1. Phase probability distribution of right circular and left circular polarization for a coherent state which overlap. 
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FIG. 2. Polarization interference because of right circular and left circular polarization for a coherent state. 



B. Phase probability density for a Squeezed state 



The properties of the squeezed states are dependent on the phase so we studied their phase properties. Taking 
account the polarization, the number states decomposition of squeezed states is written as 
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Where £ = re 10 squeeze parameter with — oo < r < oo and < 9 < 2ir and 7 = a cosh r + a*e tf> sinhr. 

Putting these values in phase probability density we get for the right circular polarization and the left circular 

polarization states as 
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We plot phase probability density for all the terms with phase varying from — n to +7r, taking the parameters as r = 1 
and 8 = and a = 1, we get the expressions as 
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The plots are, 

probability density(p(0)) 





A 0-5 








I \ 0-4 








0.2 








0.1 







-3-2-1 12 3 

FIG. 3. Phase probability distribution of right circular and left circular polarization for a squeezed state which overlap. 
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FIG. 4. Polarization interference because of right circular and left circular polarization for a squeezed state. 



C. Phase probability density for a thermal states 

Thermal states are defined by average photon number and have a completely random phase distributed over the 
period. The thermal states are a mixture of all the states, therefore, a thermal state is written as 

\^>=T, p n\n>=Ea^\ n >+EafL^ l - n - 1> ^ m 

n=0 y ' n=0 V ; 



where P n — (T+=b+i i is Bose- Einstein weight summed over all the states. 
Putting these values in phase probability density, we get for circular polarization and opposite circular polarization 

as, 
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and for polarization interference 
We plot the phase probability density of all the terms with phase varying from — n to +tt, taking the parameter n — 1 
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FIG. 5. Phase probability distribution of right circular and left circular polarization for thermal state, which overlap. 
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FIG. 6. Polarization interference because of right circular and left circular polarizations for the thermal state. 



D. Conclusion 



We have used the two component form of the Maxwell equations first derived by Wigner pJJ and subsequently 
studied in detail by Bialynicki-Birula, Sipe and others to construct a phase operator for photons, which is the canonical 
conjugate of the number operator. To do so we introduced an operator corresponding to the physical process of 
reversing the sense of circular polarization. Our method is sufficiently general to be extended to the case of linear 
isotropic and anisotropic media. As the states used here correspond to the two polarization states studied by Wigner 
and Kim [32] , consequently, all the operators derived by them to represent phase plates, attenuators, polarizers etc., 
can be directly applied to our case without any changes whatsoever even in the limit of vanishing photon numbers. 
It is interesting to note that it is in this limit of the double vacuum state that the Heisenberg algebra is replaced by 
complete commutation. Garrison [5 had suggested that these operators should not satisfy the commutation relations 
in the strong (Heisenberg) sense but only in the weaker (Weyl) limit as otherwise the number phase algebra will be the 
same as the position momentum algebra and this would contradict the requirement that phase has a 2ir periodicity. 
These methods are best suited for the study of phase and polarization effects in anisotropic and conducting media 
as these absorptive in nature and the elegance of the conventional Schwinger representation has to be augmented by 
phenomenological decay terms, while in the wave function representations advocated here, single boson annihilation 
operators are just as natural as the two boson operators. The other great advantage of the Bialynicki-Birula, Sipe 
representation is that they permit a position representation for photons, which is consistent with the Newton- Wigner 
theorem [39]. This permits the investigation of phenomena like EPR and other correlations [27] in a intuitive and 
physically understandable model. As such experiments are always measurements of polarization coupled with position 
our methods should find applications there too. 
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